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WEAK MATRIX ELEMENTS ON THE LATTICE: 
RECENT DEVELOPMENTS IN A'-PHYSICS 
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INTRODUCTION 

In this talk we present some recent developments in the calculation of weak matrix elements on 
the lattice. Lattice QCD is one of the few systematically improvable methods for computing 
them from first principles, and has proven a powerful and appealing approach. In spite 
of the successes, progress has been slow due to the presence of systematic effects, such as 
discretization and non-perturbative renormalization effects. In the following, we concentrate 
on the applications in A-physics of a recently introduced method for non-perturbative (NP) 
renormalization Q. 

Renormalization of lattice operators is a crucial ingredient in the calculation of physical 
weak matrix elements on the lattice. A physical amplitude A a ^,n of a weak transition a — > f3 
is calculated via the Operator Product Expansion (OPE) by 

A a ^ = C w (n/M w ){a\d(ti)\(3) (1) 

where Cw is the Wilson coefficient of the OPE, My/ is the mass of the W boson, \i is the 
renormalization scale and (a\0(fj,)\/3) is the matrix element of the renormalized operator (at 
the scale fj,) relevant to the physical process. The Wilson coefficient C\y([j,/Mw) contains 
the short-distance information and can be calculated in Perturbation Theory (PT) in the 
continuum at the renormalization scale \i. The matrix element contains the long-distance 
dynamics and thus must be calculated non-perturbatively on the lattice. Renormalization 
relates the regularized lattice matrix elements to its continuum counterpart. 

On the lattice, chiral symmetry is explicitly broken with Wilson-like fermions. The possi- 
bility of recovering the chiral symmetry in the continuum limit was shown in . The general 
prescription is to subtract from the bare operator 0{a) all the operators of dimension less or 
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equal than O(a), which have the same quantum numbers left unbroken by the regularization: 
(a\6(jj)\P) = lim(a|Z (//a)[0(a) + £ ^0,(a)]|/3), (2) 

i 

If the subtracted operators 0{ have lower dimension than O, the mixing constants are power- 
divergent in the cutoff, Z{ ~ l/a d with d > 0. These divergent factors can pick up expo- 
nentially small contributions in the strong coupling, yielding a finite contribution as a — > 0, 
i.e. 

l e -l/a.(«) „ AqcD- (3) 

a 

These divergences must be subtracted in a completely non-perturbative way. Recently, non- 
perturbatively renormalization has witnessed a great progress ||, Here, we want to 
discuss the applications of the method of ref. [[3J , which in the following will be refered to as 
the NP method (NPM). 

We will concentrate here on two physical processes of phenomenological interest: 

• K° — K° mixing, in which the bare operator AS=2 = sj^dsj^d mixes only with 
operators of the same dimension || 

From the matrix element of the AS = 2 operator the kaon 5-parameter Bk is obtained, 
which is a quantity of great phenomenological interest, being related to the e-parameter 
which measures CP-violation in the K°-K° system. With the measured value of the 
top quark mass, an accurate prediction of Bk enables us to limit the range of values of 
the CP-violation phase 5. 

• K — > 7T7T decay, in which the bare operators 0^ L also mix with operators of lower 
dimensionality 0. 

These decays are relevant to the study of AI = 1/2 rule, the enhancement of the 
AI = 1/2 with respect to the AI = 3/2 rate, a long-standing question which lattice 
QCD should be able to understand in a quantitative fashion. Yet, this has proven a 
formidable task due to the power-divergent mixing. The NPM allows us to make the 
subtraction without loosing predictive power, as one would do imposing renormalization 
conditions directly on hadronic states §|. 

NON-PERTURBATIVE METHOD 

In the NPM, the renormalization conditions are applied directly to the Green functions of 
quarks and gluons, in a fixed gauge, with given off-shell external states of large virtualities 
fim. The method mimicks what is usually done in the perturbative calculation, but the Green 
functions are evaluated in a non-perturbative fashion from Monte Carlo simulations. 

To give the flavour of the method, let us consider the simplified case of a multiplicatively 
renormalizable operator, e.g. a two-quark operator O = qTq. Given the bare lattice operator 
latt (a), the renormalization condition we impose is |Q] 

Z^^a)( P \0^(a)\p)\ p2= ^ = (p|0 latt (a)|p)|^ 2 , (4) 

where (p\ ■ ■ ■ \p) denotes the matrix element of external quarks of momenta p which can be 
calculated to all orders in the QCD coupling via Monte Carlo simulations. The renormalized 
operator obtained with the NPM is then 

d RI ( / u) = z^ I tt (Ho latt («), (5) 
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which depends on the external states and the gauge, but not on method used to regulate 
the ultra-violet divergences. To stress this point, we call the NP renormalization scheme 
Regularization Independent (RI) ||. The physical operator 

0^ B {M W ) = C m (M w /iJ,)dja(») (6) 

is independent of external momenta and gauge (up to higher orders in continuum PT and 
lattice systematic effects) if the Wilson coefficient function C-Ri(M\y / n) in the RI scheme 
is calculated with the same external momenta and gauge of Oja(fi). The advantage of the 
RI scheme is that it completely avoids the use of lattice PT, which is expected to have a 
worse convergence than the continuum expansion fpTof] . The coefficient function Cri(M\y / fi) 
are instead calculated in continuum PT, which cannot be avoided since the Wilson OPE is 
defined perturbatively. 

The NPM is valid for any composite operator, as long as we can can find a window in the 
range of fi such that Aqcd "C ^ <C 0(1 /a), in order to keep under control both the higher- 
order effects in the (continuum) perturbative calculation of Crj and discretization errors Q. 
We stress that this requirement is common to all NP methods on the lattice. 



AS = 2 

We consider the renormalization of the four-fermion operator 

O as = 2 = (srfd)(srfd) , = ~ 7m (1 - 7s) , (7) 

which appears in the weak effective Hamiltonian relevant for K°-K° mixing. 

The Wilson term in the quark action induces the mixing of the operator (0) with dimension- 
six operators of different chirality, determined by CPS symmetry of the action (the S stands 
for the s <-> d flavour exchange symmetry) |L4|, M, |j 

°' s -ik l0ss - 0pp] 

(N - 1) W 

°3 = \*at Pss + Opp + Ott] 
loW, 



O4 = ( ^g Ar 1) ' lOss + Opp - \0 TT \ 



where A^ c denotes the number of colours and Orr = (sTd)(sTd), with T one of the Dirac 
matrices. Note that O4 is not present in 1-loop PT JllJ , so the mixing is at most of C?(<7g(a)), 
where qq (a) is the bare lattice QCD coupling. We define the renormalized operator as 

Ok? =2 (p) = Z* s =\»a)0£t\H ee z^= 2 (0 A5 = 2 + j: ZiOi) (9) 

i=i 

Since there are no AS = 2 operators of dimension lower than six, the mixing constants Z{ are 
finite, while the logarithmically divergent Z AS=2 renormalizes multiplicatively the subtracted 
operator 0^ =2 . 

We determine the mixing constants Zi with a projection method on the four-point ampu- 
tated Green functions 0, 0]. Denoting by Oq the operator AS=2 , we define a set of mutually 
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Figure 1: Mixing constants Zi(i = 1, . . . , 4) at k = 0.1432, for several renormalization scales 
[j, 2 a 2 . The solid line is the result from "standard" PT, while the dashed line comes from 
"boosted" PT with a v ~ 1.68 a\ att . 
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orthogonal projectors Pj i = 0, . . . , 4 on the amputated tree- level four-quark Green functions 
A. ^ of the operators Oi, by the relation Tr PjA^ ' 1 = <5y, where the trace over colour and 
spin is understood. The mixing constants Z% are fixed by the condition that the subtracted 
operator 0^ =2 ^ e proportional to the bare free operator, i.e. 

Tr P fc A s A i= 2 (pa) = 0, fe = 1, ... ,4, (10) 

where A^ =2 (pa) is the amputated Green function of 0^b =2 > which is calculated at equal 
external momenta p and in the Landau gauge in a completely NP fashion from numerical 
simulations. Eq. ( |io|) yields a linear non-homogeneous system, from which we determine the 
mixing constants Z{. Once we have determined the mixing constants, the overall renormal- 
ization constant Z AS=2 is determined by ||, [J 

Z AS = 2 ^a)Z q 2 (fia)T^t 2 (pa)\ p 2 = ^ = 1, (11) 

where T^^ =2 (pa) = Tr PoA^ =2 (pa), and Z q is the light-quark renormalization constant, 
determined from the conserved vector current in a NP way Q. 

The quark Green functions have been calculated using Monte Carlo simulations. We have 
used an ensemble of 100 configurations, on a 16 3 x 32 lattice, at j3 = 6.0, and with three values 
of the hopping parameter (which is related to the quark mass) k = 0.1425,0.1432,0.1440 for 
the quark propagator, in the lattice Landau gauge. To reduce lattice artefacts, we have used 
an 0(a) improved lattice quark action [12|. We have performed the calculation for a wide 



range of scales [i 2 a 2 . In fig. |], we show the NP mixing constants at K = 0.1432 as a function of 
/x 2 a 2 , and compare them with the PT result in the same gauge and external momenta [Bj. We 



also report the PT result using a "boosted" coupling ay ~ 1.68 a^ att [ 10 1 . We have checked 
that the dependence of the Z's on n is very mild. As expected the higher-order contributions 
differentiate the mixing constants with respect to the PT value, which at one-loop is the 



same for all the constants [11]. We note that Z2 and Z\ are very well defined and almost 
scale independent in a large "window" of (i 2 a 2 , whereas Z\ and Z% are more scale-dependent. 
Moreover, Z^ which is absent in 1-loop PT is not neglegible. 

The effects of the NP corrections can be most clearly seen in the study of the chiral 
behaviour of the matrix element {K®\0 AS=2 \K®)\ &tt . The bare matrix elements have been 
computed from an ensemble of 460 configurations with the same 0(a) improved action, on a 



18 3 x 64 lattice, at (3 = 6.0 [13]. Parametrizing the lattice matrix element as 

(K°(p)\d&?= 2 \K°(q)) = a + 0m 2 K + 7 (p • q) + ... (12) 

one finds the values of a, (3 and 7 at several renormalization scales /i 2 a 2 , cf. tab. |l[ The 
Bk parameter in the RI scheme is given by = j(fjL)/Z^, where Za = 1.06(3) is the axial 
current renormalization constant, a is a lattice artefact and should vanish in the continuum 
limit. 

As can be seen, in the window of fi 2 a 2 values in which the Z's were found to be reliable, 
a is compatible with zero within one standard deviation, whereas (3 is compatible with zero 
within 1.5 standard deviations. On the other hand, in the same window of n 2 a 2 , 7 stabilizes 
to a non-zero value. The fluctuations of Z\ and Z3 do not influence the stability of the 
results since the relative matrix elements (K°(p)\Oi^\K°(q)) are quite small compared to 
(K (p)\02\K (q)), whose mixing constant Z2 is extremely well determined. The vanishing of 
a allows us to conclude that the use of the NP Z's improves the chiral behaviour for a large 
range of values of fi 2 a 2 ||. 
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fi 2 a 2 


a 


P 


7 


0.31 


0.030(18) 


0.27(21) 


0.90(15) 


0.62 


-0.027(16) 


0.36(18) 


0.75(13) 


0.96 


-0.012(14) 


0.24(17) 


0.69(12) 


1.27 


0.005(13) 


0.14(16) 


0.68(12) 


1.39 


-0.009(13) 


0.24(16) 


0.67(12) 


1.85 


-0.003(13) 


0.18(16) 


0.66(11) 


2.46 


-0.001(12) 


0.24(15) 


0.65(11) 


4.01 


-0.002(12) 


0.44(15) 


0.67(11) 


BPT 


-0.052(12) 


0.16(15) 


0.62(11) 



Table 1: Values of a,P and 7 for several renormalization scales /U 2 a 2 . In the last row we 
present results obtained with "boosted" PT renormalization constants. 



One we are confident that the correct chiral behaviour is recovered, we can reliably 
calculate the physical value of Bk , obtained by calculating the Wilson coefficient in the same 
RI scheme g3|,§. 



AI = 1/2 

In the continuum, with an active charm quark and the GIM mechanism at work, the operator 
basis is given by || 

°i L = \[{sd) L {uu) L ± isu) L {ud) L ] - (u - c). (13) 

Again, the renormalization strategy is complicated by chiral symmetry breaking. In fact, 
the Wilson term induces the mixing of O ll with lower-dimensional operators, with power- 
divergent coefficients, which need to be subtracted non-perturbatively. In order to renormalize 
the operators (|i~3|) on the lattice we find it convenient to separate the {8, 1} and the {27, 1} 
components under the SU(3)l <8> SU(3)# chiral group. In the following, we shall concentrate 
only on the octet component of O ll , which we will denote with Oq |?]][]. 

There are different ways of calculating the K — * tttt matrix elements, which correspond 
to different renormalization structures. We considering here a general structure of the form 



o ± = Z 



± 



o^ + Y^zfot + zfo. + zfo. 



(14) 



refering to JjJ for a more detailed analysis. In eq. (14) are the bare operators, Of : , i 



'0 

1, . . . ,4 are dimension-six operators of wrong chirality (similarly to the AS = 2 case), O5 
is a dimension-five of the form sYj ^yF^yd (S^j, = 07^ or a^v) and O3 is a dimension-three 
operator of the form sTd (T = 1 or 75). 

According to the NPM, the mixing Z's are determined by finding a set of projectors 
on the tree-level amputated Green functions (GF), with off-shell quark and gluon external 
states, the choice of which depends on the nature of the operators at hand. For the AI = 1/2 
operators we choose the following set of external states: qq, qqg, qqqq, with the momenta given 
below in eq. (|l5|). For each choice of external states, i.e. for each different set of GF, we need 

1 The lattice penguin operators, being proportional to (mj: — m^)a 2 <s!C 1, will be neglected in the following. 
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different type of projectors. Let us denote with P 3 the projector on the qq GF of the operator 
O3, with P5 the projector on the qqg GF of the operator O5, and with Pj, j = 1,...,4 
the set of mutually orthogonal projectors on the operators Oj, i = 1,... ,4. We refer the 
reader to refs. [fj], ^] for the explicit expressions of the projectors. Applying the projectors to 
the corresponding NP GF of the renormalized operators O^, with an appropriate choice of 
the external states, we require that the renormalized operators be proportional to the bare 
operators, O^Qu) oc O (a) (up to terms of 0(a)), i.e. we impose the following renormalization 
conditions (trace over colour and spin is understood in the projection operation): 

P 3 («(p)|O ± |g(p)) = 

P 5 (q(p-k)g(k)\O ± \q(p))=0 (15) 

pfiqipMp^d^qipMp)) = 0, j = 1, ... ,4 

where p and k denote the momentum of the external quark and gluon legs. The system 
of equations (|l5|) completely determines in a NP way the renormalization constants, as we 
have six conditions (non-homogeneous due to the matrix elements of Oq, cf. eq. (|i~4|)) in six 
unknown mixing constants, i = 1, . . . , 4, Z^ 1 , Z 3 . 

Unfortunately, since solving eq. (|l5|) involves delicate cancellations between large contri- 
butions, it may very likely result in a very noisy determination, even with large statistics. An 
equivalent strategy we can adopt is: 

1. We introduce an intermediate subtraction for the dimension-five and -six operators 



Of = Of + cf' 4) 3 , i = 0,...,4, 
3 , 

(±.0 ,^ 



o 5 = o, + ci 5) o 3 , 



(16) 



and determine the power-divergent mixing constants C3 ' and C3 by imposing 

P 3 (q(p)\Ot\q(p)) = 0, i = 0,...,4, 

P 3 (g(p)|o 5 |g(p)> = o. 



(17) 



2. The finite mixing constants Z i and Z 5 , which in principle can be calculated in perturbation 
theory (PT), are then determined from the system 



P 5 (g(p - k)g(k)\ Ei ZfOf + Zf6s\q{p)) = P 5 (g(p - k)g{k)\Of\q{p)) 
Pj(Q(p)q(p)\ Ei ZfOf + Z±0 5 \q(p)q(p)) = ■p j (q(p)q(p)\0±\q(p)q(p)), j = 1, . . . ,4. 



The numerical implementation of this program is in progress 



CONCLUSIONS 



Over the last few years there has been considerable progress in the application of lattice 
QCD to calculate weak matrix elements. Much control has been gained over the two major 
systematic errors one had to face, discretization effects and higher orders effects in the lattice 
perturbative expansion. Discretization effects have been greatly reduced by the "improve- 
ment" program [l!| , while the results obtained with the non-perturbative renormalization 
program make us confident that we can realiably calculate weak matrix elements in the low- 
energy regime. 
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